function res = check_parity( vec )
    res = mod( sum(vec), 2 );
end
